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Abstract 

Quantum phenomena of interest in connection with applications to computation and commu- 
nication almost always involve generating specific transfers between eigenstates, and their linear 
superpositions. For some quantum systems, such as spin systems, the quantum evolution equation 
(the Schrodinger equation) is finite-dimensional and old results on controllability of systems defined 
on on Lie groups and quotient spaces provide most of what is needed insofar as controllability of 
non-dissipative systems is concerned. However, in an infinite-dimensional setting, controlling the 
evolution of quantum systems often presents difficulties, both conceptual and technical. In this 
paper we present a systematic approach to a class of such problems for which it is possible to 
avoid some of the technical issues. In particular, we analyze controllability for infinite-dimensional 
bilinear systems under assumptions that make controllability possible using trajectories lying in a 
nested family of pre-defined subspaces. This result, which we call the Finite Controllability Theo- 
rem, provides a set of sufficient conditions for controllability in an infinite-dimensional setting. We 
consider specific physical systems that are of interest for quantum computing, and provide insights 
into the types of quantum operations (gates) that may be developed. 



I. INTRODUCTION 



Over the last three decades, there has been a steady stream of papers in the 
physics/chemistry literature [U El El lU E] that describe new experiments in atomic and 
molecular science, and new ways of thinking in which control-theoretic ideas are of central 
importance. More recently, the driving force has been the desire to manipulate quantum 
states in ways that would make possible quantum computation or quantum communication, 
see for example [61 El El Ell EH]- Phenomena involving the interaction between electro- 
magnetic radiation (light) and matter (e.g. ions, spin states, etc.) are especially interesting 
because they are possible paradigms of future quantum computing devices [12]. Many of 
the exciting ideas are related to the control of these systems (see for example, work on 
control of trapped- ion quantum states [131 US H21 EEl El HE]). Some infinite-dimensional 
systems can be made to be effectively finite-dimensional by either bandwidth limits imposed 
by the control fields [6] , or by turning off specific transitions in order to truncate the Hilbert 
space [T7j , and the controllability of such systems can be analyzed using finite-dimensional 
methods PH EDI ED E2] • 

We are interested in the quantum systems that are modelled as finite-dimensional for 
quantum computing purposes, when in fact they are infinite-dimensional. The well-known 
paper by Huang, Tarn and Clark [lj seems to be overly pessimistic with repect to the 
control of infinite quantum systems by asserting that "using piecewise-constant controls, 
global controllability cannot be achieved with a finite number of operations". In 2000, 
Zuazua summarized the field aptly thus: "From a mathematical point of view this [when 
the finite-dimensional system approaches the infinite-dimensional Schrodinger equation] is a 
very challenging (and very likely difficult) open problem in this area" [23J. Recently, Turinici 
and Rabitz [21] adapted the classical Ball-Mar sden-Slemrod results [25] to the quantum 
setting and showed that exact controllability does not hold in infinite-dimensional quantum 
systems (see also [2B]). Along with an excellent review, Illner, Lange and Teismann [2T] 
show that the Hartree equation that is well-known in quantum chemistry with bilinear 
control is not controllable in finite or infinite time. In spite of these negative results, most 
quantum computing systems that are indeed infinite-dimensional have shown themselves to 
be remarkably amenable to the production of a variety of states. Is it possible then to make 
a statement regarding the reachable set of states in such systems? 
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There has been recent progress in attacking this problem, and a few positive results. 
Ref. [2H] showed that for an infinite-dimensional quantum control system with bounded 
control operators a quantum state could be steered only within a dense subspace of the rele- 
vant Hilbert space. In 2003 [16], we showed that one can reach any finite linear superposition 
of states in a two-level system coupled to a harmonic oscillator by the alternate application of 
control fields even when one of the control operators is unbounded. Ref. [29] presents an in- 
genious method of creating arbitrary finite superpositions of rotational states of molecules by 
ensuring the trajectories lie within SU(N)[z,z-l] loop group of N-periodic matrices. Ref. [3U] 
presents a scheme that is well-known in atomic physics - when an infinite-dimensional quan- 
tum system has unequally spaced bound state energy levels, single resonant fields can be 
used to transfer population (albeit very slowly). Ref. [31] provides a prescription for deter- 
mining whether there exists a submanifold that is strongly analytic controllable, but this 
prescription does not identify the submanifolds. Ref. [32] presents an adiabatic method of 
controlling a sequentially connected system in which the transition couplings get weaker as 
one moves away from the ground state. Ref. |33j presents an algebraic framework to deter- 
mine the nature of controllability of some infinite-dimensional quantum systems, specifically 
those with continuous spectra. Ref. [31] proves approximate controllability for bound states 
of quantum systems that are not equally spaced in energy even when the control matrices 
are unbounded. 

We begin this paper in Section [IT] by discussing the difficulties associated with applying 
traditional methods of analysis to infinite-dimensional control systems. We also discuss a 
specific class of infinite-dimensional systems — one in which it is desirable to steer a fi- 
nite superposition of eigenstates to another finite superposition of eigenstates. We present 
a general theorem on Finite Controllability of quantum systems that exhibits the condi- 
tions needed for such transfers. For the specific quantum-computing models of interest, 



the trajectories are constrained to lie within subspaces. Section III presents examples of 
quantum-computing systems such as a model of a trapped-ion qubit and trapped-electron 
qubits. We show how our theorem can used to be prove finite controllability in the former 
case. We also discuss other similar systems where the theorem can and cannot be applied. 
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II. INFINITE-DIMENSIONAL CONTROLLABILITY 



Controllability results for infinite-dimensional systems are seldom just straightforward 
extensions of the finite-dimensional ones, and in particular this is true for bilinear sys- 
tems. Recently, there has been significant interest in the class of bilinear systems because 
of their relevance to quantum control. In the following subsection we illustrate the lim- 
itations of applying the tools of finite-dimensional systems analysis to certain classes of 
infinite-dimensional systems. 



A. Limitations of Lie Algebraic analysis 

Examining the Lie algebraic structure often gives us insights into the controllability of a 
quantum system, but in the case of infinite-dimensional systems, this insight is limited. In 
the well-known example of a resonantly-driven quantum harmonic oscillator (see e.g. [35J, 
[T6] and [36], the evolution is given by 

(&-■>)-*«•)* (1) 

Here, the bilinear control term u(t)x arises because of the dipole interaction between the 
field and harmonic oscillator. The two operators of interest, A = | — x 2 ^j and B = —ix 
generate a Lie algebra of skew-hermitian operators that is just four- dimensional. Thus, we 
expect that the control of this system will be limited. For example, it is well-known that it 
is not possible to transfer the number state x(0) = |0) to x(T) = \n) for n > [37] using 



this control. We return to this example in Section [Hi 

However, even when we encounter infinite-dimensional systems for which the Lie algebra 
also is infinite-dimensional (such as in Ref . [TBI [38] ) . the statements one can make about the 
controllability are also limited. More work is required to say with precision exactly what 
the reachable states are. This issue is illustrated in Section IIII1 

In computing the span of the Lie algebra, it is necessary, of course that the domain of 
the operators involved be such that the bracket operations are allowable. This is the case 
in this paper, but we do not discuss these technical details here. 
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B. Finite Controllability 



In this section, we prove an elementary but useful theorem about controllability on finite- 
dimensional subspaces of a complex Hilbert space. It is the basis for the applications we 
will present in the next section. 

Definition 1 We will say that a finite- dimensional system evolving in the space of complex 
n-vectors x, 

x = UiGiX 

i 

with skew-Hermitian operators G{, is unit vector controllable if any unit length vector xo 
can be steered to any second unit length vector Xf in finite time. 

Theorem 2 (Finitely Controllable Infinite Dimensional Systems) Consider a com- 
plex Hilbert space X together with a nested set of finite- dimensional subpaces H = {Hi C 
7i 2 C Hz-}- Consider 

m 

% = C^UiG^x. 
i=i 

Assume that Hi is an invariant subspace for a subset Q\ of the set {Gi} and that the system 
is unit vector controllable on Hi using only this subset of the Gi. If for each H a 1 there 
is a subset Q a of {Gi} that leaves H a invariant and if for any unit vector in H a the orbit 
generated by exp (Q a ) contains a point in one of the lower dimensional subspaces Hp then 
any unit vector in any of the Hi can be steered to any other unit vector in any other Hj 
using a finite number of piecewise constant controls. 

Remark: Given a system and a nested set of finite dimensional subspaces it will be said to 
be finitely controllable if it can be transfered from any point in one of the subpaces to any 
other point in that subspace with a trajectory lying entirely within the subspace. 

Proof: Suppose that x(0) G H a and x(T) G Hp are unit vectors representing the initial 
value of x and the desired final value, respectively. Then by assumption either there exists a 
subset of the {Gi} that leaves H a invariant and steers x(0) to a point in some Hp with the 
dimension of Hp being strictly less than that of H a or else x(0) G Hi. A finite induction 
on the index of the set Hi then shows that the state can be steered to a unit vector in the 
subspace Hi, which is a controllable subspace. To finish the proof, observe that if one can 
reach x a G H from x fe G H then the standard time reversal argument using the fact that 
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the Gi are skew-Hermetian shows that it is possible to reach Xb G Ti from x a G Ti- Thus a 
second application of the steps given above implies that it is possible to reach x(T) G Tip. 

Discussion and example: Various technical issues that arise in infinite-dimensions illus- 
trate the significance of this theorem. All finite-dimensional subspaces are closed in the 
usual topology but infinite-dimensional subspaces need not be and there are a number of 
issues involving infinite-dimensional controllability that can be clarified by giving more at- 
tention to the distinction. Let I2 denote the Hilbert space of infinite vectors whose entries 
are square summable which corresponds to the entire state space of our quantum system. 
Let l denote the subspace consisting of those elements with only a finite number of nonzero 
entries, which corresponds in our setting to a finite superposition of states. Of course, this 
subspace is not closed. 

We recall now some general ideas from operator theory [39] and discuss their relationship 
to our theorem. 

1. If e and e B are bounded operators mapping a Hilbert space Ji into itself and if V is 
an invariant subspace for e A and e B , then V is invariant for the product e A e B . 

In our setting e Bc and e Br are unitary and thus bounded. 

2. If A is the infinitesimal generator of a semigroup mapping a Hilbert space Ji into itself 
and if V is a closed, nontrivial, invariant subspace for A, then V is invariant subspace 
for e At . 

Consider for example an operator A on I2 that maps the i th unit basis vector e; into 
ej + i for all i = 1, 2, .... This operator clearly sends Iq into itself but e At for t 7^ sends 
ei into the element Z^e^/i! which is not in l . 

3. If A and B are infinitesimal generators of a semigroup mapping a Hilbert space Ti into 
itself and if V is a closed, invariant subspace for A and B, then V is invariant subspace 
for e ( A+B )'. 

Again, the assumption that V is closed cannot be dispensed with. For example, con- 
sider a system x = uAx + vBx with A being block diagonal with two-by-two blocks 
and B being basically block diagonal with two-by-two blocks except for top and left 
borders where the necessary one- by-one, two-by one and one-by two matrices pad an 
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otherwise block diagonal matrix. Specifically, A and B are given by 



.4 



B 



Bq O22 O22 O22 •• 

O22 Bq O22 O22 •• 

O22 O22 -Bo O22 •• 
O22 O22 O22 B .. 

On 12 12 0i 2 ... 

O21 B O22 2 2 ••• 

O21 O22 B O22 ■•• 

O21 O22 O22 Bq ... 



with Oij denoting an i by j matrix of zeros and 



Bn 



1 
-1 



(2) 



(3) 



(4) 



The operators A and B leave Iq invariant but the exponential of their sum does not. 

One can provide a direct argument not involving Lie theoretic techniques that shows that 
any unit vector can be transferred to e\. The idea, motivated by the analysis in Ref. [T3] , 
is to alternate the use of A and B and reason that if we start with a vector with x n 7^ 
and x m>n = we can use a control with v — (or u = depending on whether n is even 
or odd) to reduce the vector to one for which x n is zero, then use a control with u = to 
reduce the vector to one with x n _i = without changing x n , etc. This argument depends 
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on observing that 



e B u 


o 22 o 22 


o 22 


o 22 


e B ° u 22 


o 22 


o 22 


22 e B ° u 


o 22 


o 22 


o 22 o 22 


e B Q u 



(5) 





1 


O12 


0l 2 


0l2 




O21 


e Bov 


o 22 


o 22 


e Bv = 


O21 


o 22 


e Bov 


o 22 




O21 


o 22 


022 





from which we see that e Au simply rotates the entries in x pairwise, i.e., Xi,a; 2 ; 23,0:4; etc. 
so that by attending only to the subspace defined by x n and one can make x n = 0. 
Having done so, we can use e Bv to rotate in x n _i,x n _ 2 space so as to reduce x n -i = with 
out changing a? m > n -i> etc. This shows that the Lie group generated by e Au and e Bv acts 
transitively on the unit vectors and thus identifies the Lie algebra to within one of a few 
possibilities. (See Ref. [20J). 

III. PHYSICAL SYSTEMS 

In this section, we apply the Finite Controllability Theorem to determine the reachable 
set of states of four infinite-dimensional quantum systems, namely, the quantum harmonic 
oscillator, a spin-half system coupled to a harmonic oscillator (model of a trapped-ion qubit), 
an N-level atom coupled to harmonic oscillator, and a spin-half system coupled to two 
harmonic oscillators (model of a trapped electron qubit). Note that in all sections below 
we use atomic/scaled/dimensionless units so that the evolution equations do not explicitly 
contain Planck's constant H, the charge of the electron e or the mass of the electron m e . 

A. System 1: Quantum harmonic oscillator 

We discuss this well-known system first. The controllability algebra is finite-dimensional 
and, in particular, the system does not satisfy the conditions needed for the application of 
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the Finite Rank Controllability Theorem. The discussion, however, is key to setting up the 
formalism used in subsequent examples that are finitely controllable. 

The problem of controlling the harmonic oscillator has been discussed many times (see 
e.g. [10] and other references as discussed above). If the control is a sinusoidal resonant 
driving field (of frequency equal to the harmonic oscillator frequency u m ) as shown in the 
transfer graph Fig. [TJ then the evolution is via 

Here, the control term u(t)x arises because of the dipole interaction between the field and 
harmonic oscillator. The operators of interest are A — ~ (Jj^ — an d B = —ix. A 
and B generate a Lie algebra of skew-hermitian operators that is just four-dimensional 
(C = [A, B] = J^, D = [B, C] = il, where / is the identity operator). This in itself tells us 
that the resonantly driven harmonic oscillator is not controllable. 



|n=o) |n=l) |n=2) | n =3) 



In=4 



FIG. 1: Graphical representation of the quantum harmonic oscillator driven by a sinusoidal res- 
onant field. Note that while the strengths of the transition couplings increase as the square root 
of the quantum number n as shown by the boldness of the connections between energy levels, the 
transition frequency between each level is the same. It is this latter property that leads to a lack 
of controllability. 

As is well-known, the spectrum of A is discrete. If we describe the evolution in terms of 
an eigenfunction expansion, with the basis being |n)'s, the eigenfunctions of d 2 /dx 2 — x 2 , 
then the evolution is via 

x n = -iu m {n+ ^)x n (8) 
-u(t)-j=(y/n - - \/nx n+1 ). 

Although the eigenstates of the harmonic oscillator can be written as an infinite set of nested 
finite subspaces, it is seen that the operator B connects space Hi to both Tti-i and TCi+\. 
Thus finite superpositions of eigenstates may not be reached by resonantly driving the har- 
monic oscillator, consistent with the fact that the requirements of the Finite Controllability 



Theorem are not met. Physically, this is due to the degeneracy of spacings between the 
eigenstates and the fact that the control vector field simultaneously illuminates all states. 

However, the system is in fact partially controllable. One can be quite explicit about the 
reachable wave functions under the application of a control [37J. It is interesting that the 
finite-dimensionality of the algebra that makes the system uncontrollable is also the property 
that makes it possible to derive an analytic expression for the unitary propagator jH]. For 
example, if we agree to call any function of t and x which has the form 



i-gaussian , then if the initial value of ip is i-gaussian the solution will be i-gaussian for all 
time, regardless of the choice of u(t). (This is directly analogous to the fact that the solution 
of the conditional density equation of estimation theory remains gaussian if it has a gaussian 
initial condition.) An example of such a function is the well-known coherent state (12] 



n=0 

We see that it is not possible to transfer x(0) = \a) to x(T) — \i) for i > because \i ^ 0) 
is not i-gaussian, but |a) is, and is in fact a coherent state. 

B. System 2: Spin-half particle in a quadratic potential 

In contrast to the harmonic oscillator, the model of a spin-half particle coupled to a 
harmonic oscillator with suitable controls turns out to be finitely controllable. This model 
is a good representation of an ion with two essential internal states trapped in a quadratic 
potential. We show below that this system satisfies the conditions of the Finite Controllabil- 
ity Theorem of Section |2j Moreover, one can also provide an algorithm for explicit control. 
The spin- 1 model represents a two- level atomic ion with an energy splitting huo, where the 
frequency ujq/2tt is in the several GHz range. The atomic levels are coupled to the motion 
of the ion in a harmonic trap [13] • These quantized vibrational energy levels are separated 
by a frequency u m /2ir in the MHz range. 

In a frequently cited paper, Law and Eberly [13] showed that when properly interpreted, 
this system has interesting controllability properties, quite different from the properties of 
the harmonic oscillator alone. In fact, by coupling the harmonic oscillator with a two- level 




(9) 




(10) 
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system it is possible to arrive at a system which is much more controllable than the harmonic 
oscillator. At an intuitive level, this can be seen simply as a consequence of the fact that 
the addition of a spin degree of freedom breaks the infinite degeneracy associated with the 
harmonic oscillator and allows the system to resonate with more than one frequency. This 
allows the transfer of population from any eigenstate to any other eigenstate by sequentially 
applying the two frequencies. We now analyze this system from a controllability viewpoint. 

An eigenstate of the spin-half system coupled to a quantum harmonic oscillator is denoted 
by \S,n), where the first index refers to the "spin" state of the system, and the second index 
is the number state of the harmonic oscillator. An applied field causes transitions between 
the eigenstates of the coupled spin-oscillator system. A monochromatic field of angular 
frequency u = u Q causes resonant transitions between states | l,n) and | |,n) (carrier 
or spin- flip transitions). A monochromatic field of angular frequency to = ujq — u m causes 
resonant transitions between states | j, n) and | n— 1), i.e., produces so called red sideband 
(that is with angular frequency to = u — uj m ) transitions. 

These transitions are graphically depicted in Fig. [2] with the thickness of the edges qual- 
itatively representing the strength of the coupling between the states. As pointed out in 
Ref. |18j . when both fields (carrier and red sideband) are applied simultaneously, the eigen- 
states of the system are sequentially connected. Therefore, we look at the trapped-ion model 
controlled only by these two fields. 



FIG. 2: Graphical representation of the coupled spin-half quantum harmonic oscillator system 
driven by sinusoidal resonant fields of angular frequency uo c and co r as shown. When t/< 1, the 
strengths of the to c transition couplings are independent of the harmonic oscillator quantum number 
n, whereas the strengths of the tu r transition couplings increase as the square root of n as shown by 
the boldness of the coupling lines. Note that there is no direct coupling between two consecutive 
oscillator states with fixed spin. 
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Now we write the evolution equation of the spin-half coupled to harmonic oscillator 
driven by two fields that drive the carrier and red sideband transitions. The amplitudes 
corresponding to the fields that cause the carrier and red transitions are dubbed E c and 



E r respectively. As detailed in Appendix IV, in the interaction picture and in the energy 
eigenbasis, the evolution equation is written as 



Y = {u{t)B c + v{t)B r )Y. 

The controls u{t) and v(t) are related to the applied fields via the equations 

u(t) = Cl E c (t) =0.25 fiexp(-r] 2 / 2) E c {t), 
v(t) = c 2 E r (t) = 0.25r]fiexp{-r] 2 /2)E r {t). 



(11) 



(12) 
(13) 



Here T], the so-called Lamb-Dicke parameter, is the product of k, wave vector of the light, 
and xq, the amplitude of the zero-point motion of the particle in the harmonic potential (or 
the spatial extent of the ground state harmonic oscillator wave function). By ordering the 
eigenstates as | |, 0), | f, 1), . . . , | j, 0), | |, 1), . . ., the control matrices are written as 

/ 




The upper-triangular matrices L and L\ are defined as 

( L {r?) 

o W) o 

L 2 ( V 2 ) 

V : : : 

o 4V) o . 
o o • 



V i i i 



(14) 
(15) 



(16) 



(17) 



/ 



This structure of the control Hamiltonian precludes the use of adiabatic methods such 
as those used in Ref. [32]. Unlike the situation encountered in the analysis of the quantum 
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harmonic oscillator algebra, here the formal commutation of the operators B c and B r does 
not lead to a finite-dimensional algebra, suggesting that the model with spin is much more 
controllable. This is the will be explored in the next subsection. 



1. Controllability: Lie Algebra 



It is interesting to compare our analysis with the formal calculations suggested by Lie 
theory. The first thing to do is to determine the formal structure of the Lie algebra, which 
we now consider. 

The control of the trapped-ion system is often studied in two different limiting cases - one 
in which the extent of zero-point motion of the spin-half particle in the harmonic potential 
Xq is much smaller than the wavelength of the applied light 2ir /k, i.e., r] <C 1 (the Lamb- 
Dicke limit), and the other in which 77 ~ 1 (beyond the Lamb-Dicke limit). The case in 
which 77 ~ 1 is more general than the case of the Lamb-Dicke limit, but requires a more 
sophisticated analysis. We study initially the Lamb-Dicke limit in which the Lamb-Dicke 



parameter ?| < 1. The terms in equations (16) and (17) are expanded to first order in 77. 



The control Hamiltonians can then be expressed in operator form as 



B r 



i 

1 



and 



B r =V 



a 
-a 1 " 



(18) 
(19) 



where a and 0) denote the annihilation and creation operators of the harmonic oscillator 



as defined in Appendix IV (Note that B r is the same Hamiltonian as obtained from the 
well-known Jaynes-Cummings model [2] that describes the interaction between a quantized 
cavity field and a two-level atom.) 

In order to compute the Lie algebra, let us consider T, an operator acting on a complex 
Hilbert space. We associate with T a skew-hermitian operator acting on Ti. © TC defined by 



J(T) 



T 
-Tt 



(20) 
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For convenience, let K{T) be another operator denned in a similar way as 



T 



K(T) 



(21) 







T 



Of course, K(T) is skew-hermitian if and only if T is. The control operators we are interested 
in for the purposes of determining the structure of the Lie algebra are given by B c = J(il) 
and B r = r]J(a). We have 

Lemma 3 The Lie algebra generated by J(il) and J(T) includes the operators 



where, W = i(T + T*). 

Proof: A calculation shows that [J(T),J(iI)} = K(W) and further, [J(iI),K(W)\ = 
—2iJ(W). We can then check that 



These calculations make it clear that if the powers of W are independent then J(il) and 
J(T) do not generate a finite-dimensional algebra. Thus if T is nonzero only on the diagonal 
immediately above the main diagonal (which is true for the operator a), and if every term 
on this upper-diagonal is nonzero, then the successive powers of W are independent and the 
algebra is infinite-dimensional. 

This is the case for the coupled spin-half harmonic oscillator system. Of course, this 
calculation only shows that this system, unlike the harmonic oscillator, does not generate 
a finite-dimensional controllability Lie algebra. More work is required to say with precision 
exactly what the reachable states are. This is precisely the role of Theorem [2] which gives 
more specific information of which operators play a role in the control process. 

Note: In the case where the Lamb-Dicke limit does not apply, the Lie algebra will still 
be infinite-dimensional but the terms are more complicated. 

2. Finite controllability 

In this subsection we discuss how finite controllability works in this infinite-dimensional 



J(W 2p ); p= 1,2,3,.. - ; 



K(W 2p+1 



); p = o,i,2,-.- , 



(22) 




(23) 



setting. 
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From Fig. |2j it is seen that the sequentially connected eigenstates can be looked at as 
an infinite set of finite-dimensional subspaces with the ground state | j,0) being equal to 
Further, when operators B c and B r are applied sequentially, each subspace Hi can 
be transferred to Thus the criteria for finite controllability are met. By sequential 

application of the two operators, any finite superposition of eigenstates can be transferred 
to the ground state in finite time. 

The application of these statements to the spin-half in quadratic potential example is 
best understood by writing the control matrices B c and B r in a re-ordered basis as follows: 
The eigenstates can be ordered as | |, 0), | |, 1), . . . , | j, 0), | j, 1), .... In the interaction 
picture, the Schrodinger equation is written as 



Y = (u(t)B c + v(t)B r )Y, 



where u(t) and v (t) are defined as before. Then 



B r 



L 
Lo 










...\ 






Li 
Li 
















L 2 
L 2 













(24) 



B r 



f° 
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rW 
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V '■ 












■■) 



Li's and L\ s are Laguerre polynomials of the zeroth and first order, all with argument rf . 
One can immediately see the applicability of the general statements made in the sub- 



section |II B| to this system. The model of the trapped-ion qubit highlights the existence 
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of important examples for which it is desirable for the evolution to occur on a non-closed 
subspace of a Hilbert space, i. e., the space of finitely nonzero elements. In this model, this 
non-closed subspace consists of vectors in the oscillator representation with finitely many 
nonzero elements. The B c and B r operators and their one parameter groups, leave invariant 
the subspace Iq of I2 consisting of finitely nonzero sequences. The semigroup e aBc+f3Br will 
not, however, typically have any nontrivial invariant subspace. To satisfy the requirements 
of the Finite Controllability Theorem one never uses a linear combination of the operators. 
Further, as we have seen, the key to controllability is that each operator has a different 
invariant subspace within the set of finite superpositions. 



3. Explicit finite controllability scheme 

The property that both control vector fields are never used simultaneously is exploited 
by Law and Eberly [13J and Kneer and Law [TJ] in order to devise a explicit scheme for the 
production of a finite superposition of eigenstates from another finite superposition in the 
control of a spin-half particle coupled to a harmonic oscillator (in the Lamb-Dicke limit). It 
shows that if x can be transferred to y by a series of such "single nonzero u" moves then 
the transfer from y to x is also possible. 

Specifically the Law-Eberly scheme [T3] to transfer any eigenstate \i) to any other eigen- 
state \j) involves the alternate use of transitions generated by spin reversal (7r-pulses of 
E c ) and transitions generated by 7r-pulses of E r which convert from a state in which the 
oscillator has energy Ei and spin down to a state in which the energy of the oscillator in 



altered by one unit and the spin is flipped as well (see equation (19)). For example suppose 
we wish to drive a state from the | |, n) to | f, n — 2) (see Fig. [2]). This can be done using 
B r to drive the system from | |, n) to | f, n — 1), B c to drive the system from | |,n — 1) to 
I I, n — 1) and finally B r to go from | j, n — 1) to | |, n — 2). 

We note that this scheme works both in the Lamb-Dicke limit and beyond the Lamb- 
Dicke limit. In the Law-Eberly scheme, the 7r-pulses of E c are all of the same time duration 
because in the Lamb-Dicke limit, all the carrier transitions are equally strong. However, the 
coupling strengths of the red-sideband transitions are proportional to y/n, and therefore the 
7r-pulses of E r are shorter in duration as eigenstates of higher n are addressed. In order to 
generate an arbitrary superposition of a finite number of eigenstates, starting from another 
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arbitrary superposition, an additional trick is to go through the ground state of the system 
which acts as a "pass state" ■ It is possible to provide an explicit algorithm which will 
drive the system from any finite superposition to any other finite superposition. 

To prepare an arbitrary finite superposition, the simplest path is to take the system 
through the ground state. One assumes that the desired state is the initial state and then 
designs a sequence of alternating pulses of the E c and E r fields that would take this state 
to the ground state | |,0) [14] . The actual sequence that produces the superposition is 
the time-reversed sequence that was designed. For example, if the desired superposition is 
(I T; 3) + | |, 2})/v2, the sequence of pulses that will transfer this state to the ground state 
is £ c (1) (tt) E^ 2 \(j) 2 ) E^((f) 3 ) E?\fa) Ef\fc) ^ 6) (0e) E<p{fa). The action of each pulse is 
the following: EfP is a 7r pulse of the carrier field that moves the state | |,3) to | |,3). 
(Simultaneously, the population in | |,2) is transferred to | 1,2)). Ef> is a pulse of the 
red-sideband field that moves between the states | |,3) and | |,2). Since there is already a 
superposition of the two states, the duration of the red-sideband field is shorter than that 
of a 7r-pulse. Simultaneously, a superposition of | 1,2) and | is created. The next 
transition E^\(f) 3 ) transfers population between | f, 2) and | 1,2), and again is shorter than 
a 7r pulse. This sequence progresses till all the population is in | |,0). The actual sequence 
is the time-reversed sequence of the one that is described above — this creates the desired 
superposition from the initial ground state. 

If one were to transfer an arbitrary initial superposition to an arbitrary final superposition 
of eigenstates, one employs the above algorithm twice. The sequences A and B that take 
the system from the initial and final superpositions respectively to the ground state are first 
calculated. Then the sequence A is first applied taking all the population to the ground 
state. The time time-reversed sequence of B is then applied which takes the population to 
the desired final superposition. Clearly, this scheme works in finite time only if the initial 
and final states are both superpositions of a finite number of states. 

Note that finite superpositions are dense in the Hilbert space of all possible states. Hence 
from our Lie algebra analysis and the use of the Law-Eberly algorithm we have 

Proposition 4 The span of the Lie algebra generated by the operators B c and B r for the 



quantum control system in Eq. (11) is infinite- dimensional and the reachable set, which is 



dense in the Hilbert space of all states, includes all finite superpositions. 
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Note that this provides an explicit dense subspace controllability result which is hard to 
prove by abstract methods (see [23] and [23J ) . Note also that the proof of controllability 
that Law and Eberly give of what they term "arbitrary control" might be more accurately 
described as demonstrating that any state in Iq can be mapped to any other state in Iq, 



staying within l (see the discussion in subsection II B). 



4- Red and blue sideband controlled trapped-ion qubit 

The importance of the requirement of the Finite Controllability Theorem that the the 
system is unit vector controllable on Tii is seen by examining the trapped-ion qubit driven 
by the red and blue sideband fields alone. Recall that a monochromatic field of frequency 
to — ujq — oj m (lj = u + u) m ) causes resonant transitions between states | j, n) and | |, n — 1) 
(| l,n) and | f, n+ 1)), i.e., produces red (blue) sideband transitions. The Schrodinger 
equation in the interaction picture is written as 

Y = ( Vl (t)B r + v 2 (t)B b )Y, (25) 

where vi(t) and V2(t) are of the form of v(t) defined before. Then in terms of the Laguerre 
polynomials of the first order and with argument rj 2 , 
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FIG. 3: Graphical representation of the coupled spin-half quantum harmonic oscillator system 
driven by sinusoidal resonant fields uj r and u;& as shown. The strengths of the u r and cj& transition 
couplings increase as the square root of n. 

From Fig. [3j it is seen that the sequentially connected eigenstates can be looked at as an 
infinite set of finite-dimensional subspaces with the space consisting of | j, 0) and | |, 0) being 
equal to 7i\. Further, when operators B r and Bj, are applied sequentially, each subspace TCi 
can be transferred to TCi-%. However, by sequential application of the two operators, any 
finite superposition of eigenstates cannot be transferred to the ground state in finite time. 
This is because the subspace TC\ is not transitively connected and therefore not unit vector 
controllable. 

Note that this is a special case because these two control fields split the Hilbert space into 
two unconnected subspaces (denoted in the graph by the states with the grey nodes and the 
black nodes). In each of these subspaces, the Finite Controllability Theorem holds. That 
is, one might transfer finite superpositions to other finite superpositions in one subspace or 
the other. 

C. System 3: Control of an N-level ion trapped in a quadratic potential 

We note that the above example may be extended to models of an iV-level system coupled 
to a quantum harmonic oscillator. Without loss of generality, it can be assumed that the 
energy levels in the iV-level system are not equally spaced (for example atomic levels). If 
N — 1 monochromatic, resonant fields are available to couple every pair of adjacent energy 
levels, the iV-level system itself is transitively connected. It is necessary to have one more 
control field in order to make the TCi to TCi-i transition (ladder transition). There are 
multiple control schemes that are in keeping with the spirit behind the Finite Controllability 
Theorem. Consider the specific case when N = 3 (the generalization to higher N is fairly 
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obvious). The eigenstates of the coupled system are shown graphically below in Fig. |4[a). 
The resonant fields u c \ and uj C 2 that transitively connect the eigenstates of the ion are also 
shown. 




FIG. 4: (a)Top: Graphical representation of the N-level ion coupled to quantum harmonic oscillator 
driven by sinusoidal resonant fields u c i, uj C 2 and u r as shown. The strengths of the u c \ and lo C 2 
transition couplings are independent of the harmonic oscillator quantum number n, whereas the 
strengths of the uj r transition couplings depend on n. (b) Bottom: Graphical representation of 
the N-level ion coupled to quantum harmonic oscillator driven by sinusoidal resonant fields u c x, 
uj C 2 and uj r as shown. This control scheme, although not very intuitive or elegant, is based on the 
Finite Controllability Theorem. 

Consider an additional field that accomplishes the ladder transition by connecting the 
|l,n) and the \N,n — 1) states. In this case, the eigenstates of the coupled system are 
sequentially connected, and Zo subspace control proceeds exactly as in the previous exam- 
ple. The control matrices also look very similar. As before, in the interaction picture, the 
Schrodinger equation written as 

Y = ( Ul (t)B cl + u 2 {t)B c2 + v(t)B r )Y, (26) 

where Ui^it) and v(t) are defined as before. Then qualitatively, with 'X', 'Y', 'Zj's denoting 
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non-zero, real matrix elements, 
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Another control scheme that exemplifies this theorem is an additional field that accom- 
plishes the ladder transition by connecting the \2,n) and the \N,n — 1) states, as shown in 
Fig. |4](b) . Control in this case is a little more complicated, because the fields that transi- 
tively connect the iV-level system must be correctly applied in order to bring the population 
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to the \2,n) states before applying the ladder transition. The control matrices do not look 
as elegant as those in the first case, and practically this is a weaker scheme. As N increases, 
it is clear that the number of such schemes increases. However, the sequentially connected 
scheme is the one that is easiest to implement. The Law-Eberly and Kneer-Law schemes 
can be modified to provide explicit control schemes in both these cases. 



D. System 4: Spin-half particle coupled to two harmonic oscillators 

In this section we consider another paradigm of quantum computing - the trapped- 
electron qubit, which can be well-modelled as a spin-half particle coupled to two harmonic 
oscillators. We show that this is a system that is less controllable than the spin-half particle 
coupled to one harmonic oscillator. In this system, it is not possible to make a transfer 
from any finite superposition of states to any other finite superposition. This implies in par- 
ticular that the system cannot be finitely controllable. As we illustrate below, the natural 
control sequence does not allow one to remain within any one of the nested invariant sub- 
spaces spanned by the eigenstates. In addition the nested sequence of invariant subspaces 
is infinite-dimensional. 

As detailed in Ref. [IB], the energy eigenstates of a spin-half particle coupled to two 
harmonic oscillators (called the cyclotron oscillator and axial oscillator respectively) can be 
written as \nlj) where n refers to the number state of the cyclotron oscillator, I refers to 
the number state of the axial oscillator, and j refers to the spin state (up or down). The 
system is addressed by three control fields. The spin qubit is controlled using a field of 
angular frequency uo s that connects states | n I |) and \n I f). The spin-axial transition is 
controlled by field of angular frequency u sa that connects states \n I [) and \n I + 1 t)- 
A spin-cyclotron transition is controlled by a field of angular frequency u sc that connects 
states \n I I) and \n — 1 I ]). The transfer graph in Fig|5] clearly indicates that the three 
fields transitively connect all the eigenstates of the spin-axial-cyclotron system. 

As in the case of the trapped-ion, the evolution equation is clearer in the interaction 



picture and is derived in Appendix IV Writing down the control matrices explicitly is not 
trivial, and doesn't provide significant insights into the problem. Instead, one can look at 
the structure of these matrices as below. The control matrices can be written in the basis of 
trapped-electron eigenstates denoted by \nlj), where n the cyclotron state, I the axial state, 
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FIG. 5: Transfer graph of the controlled trapped-electron. States are denoted \nlj), where n is 
the cyclotron harmonic oscillator number state, I is the axial harmonic oscillator number state, 
and j is the spin state. Lines marked u SJ uisai and uj sc indicate spin, spin-axial and spin-cyclotron 
transitions, respectively. 



and j is the spin state. Qualitatively, S, A and C indicate spin, spin-axial and spin-cyclotron 
transitions, respectively, and the prefactors indicate the relative strengths. The eigenstates 
are ordered as: \l — 0, n = 0, 1), \l = 0, n — 0, |), \l — 0, n = 1, j), \l — 0, n = 1, f), |/ = 

0, n = 2,|), |i = 0,n = 2,T>, ...,|/ = l,n = 0,I), |i = l,n = 0,T>, |/ = l,n = l,|), \l = 

1, n = l,T), |/ = l,n = 2,I), |/ = l,n = 2,T), .... 

The control matrices that describes the spin-flip and spin-axial transitions have the form 
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(30) 
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respectively. 
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Here X$ and Xa are the infinite block matrices 
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respectively. 

The control matrix describing the spin-cyclotron transition has the form 
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where Xc is an infinite block matrix of the form 
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When the three transitions are applied sequentially, it is possible to transfer population 
from subspace Hi to subspace TCi-% only under certain conditions. Trivially, if one wishes 
to transfer population within one of the harmonic oscillator states, the problem is the same 
as that of the trapped ion. Also, it is easy to see that the levels are connected in such 
a way that the system is eigenstate controllable in the sense that the population can be 
coherently transferred from any eigenstate to any other eigenstate. For example, consider 
the set of eigenstates illustrated in Figj5j The condition for eigenstate controllability is that 
the pulses of frequency u> sa , uj s and u sc must be applied sequentially, and not simultaneously. 
For example, let us say we want to transfer the 1 000) state to the | 111) state. We can do so 
by the pulse sequence: p s {^) that transfers to |001), p sa {^) that transfers to 1 010) , and then 
Psc(tt) that transfers to 1 1 1 1) . 

Superficially this system appears to be very similar to the trapped-ion system. However, 
it is seen that this system is not finitely controllable; that is, even though the system is 
eigenstate controllable, it is not possible to transfer a superposition of eigenstates to the 
ground state even with sequential applications of the field. This can be seen visually from 
the control graph in FigJSJ or by examining the control matrices as follows. 

For example, let us say we want to transfer the 1 000) + 1 1 1 1 ) state to the 1 000) state. The 
pulse sequence: p s (jt), p sa (n), p sc (7r) that transfers the |000) state to the 1 1 1 1) state also 
transfers the 1 1 1 1 ) state to the |000) state. Explicitly, p s (tc) transfers 1 1 1 1 ) to 1 1 10) , but 
also transfers 1 000) to 1 00 1) ; p sa {^) transfers 1 1 10) to 1 10 1 ) , but also transfers |001) to |010); 
and Psc{k) transfers 1 101) to 1 000) , but also transfers 1 010) to 1 1 1 1) . The pulse sequence 
that transfers population down one of the harmonic oscillator ladders, transfers population 
up the other harmonic oscillator ladder. Thus we cannot apply the Finite Controllability 
Theorem using the above basis. 

This lack of controllability can be attributed to the fact that the invariant subspaces are 
not finite-dimensional. In fact this system presents an infinite number of nested infinite 
-dimensional subspaces. 

Thus we see that eigenstate controllability is a much weaker condition than subspace con- 
trollability. Simply having vector fields that sequentially connect the eigenstates is sufficient 
for eigenstate controllability. This is in contrast to the controllability of finite-dimensional 
quantum systems where eigenstate controllability implies controllability of finite superposi- 
tions. 
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This analysis provides important input into the the development of quantum gates in the 
trapped-electron system — it is not possible to use Law-Eberly type methods to construct 
arbitrary quantum gates in this system. 

IV. SUMMARY 

Many of the novel questions that arise in laying the ground work for quantum computing 
can be thought of as questions about the controllability of Schrodinger's equation. In this 
paper we prove a Finite Controllability Theorem which is useful for analyzing certain infinite- 
dimensional quantum control problems. We discuss several related physical systems that 
are models for quantum computing. 

Among the more interesting paradigms of quantum computing is the trapped ion modeled 
as a spin-half particle coupled to a quantum harmonic oscillator. In this paper, we discuss a 
general setting for this type of problem based on infinite-dimensional differential equations 
and Lie groups acting on a Hilbert space. This allows us to explore the Lie algebraic 
approach to controllability in this setting. In particular, we show that even though the formal 
Lie algebra associated with the Jaynes-Cummings model is infinite-dimensional, explicit 
finite controllability results can be determined. We also establish generalized controllability 
criteria for improving the controllability of some important infinite-dimensional quantum 
systems. 

For the specific quantum computing system of the trapped-ion qubit, we showed that the 
Law-Eberly control scheme does not require the system to be in the Lamb-Dicke limit. We 
also showed that finite controllability cannot be achieved in the trapped-electron quantum 
computing paradigm, thus limiting the type of quantum operations (gates) that can be 
executed in this system. 
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Derivation of the evolution equation for the trapped-ion qubit system 

The trapped-ion qubit is modeled as a spin-half particle coupled to a quantized harmonic 
oscillator. Since various aspects of this model are discussed in several publications, we 
summarize only the essential features here. Consider the description of a particle with 
two spin states in a quadratic potential field. The Hamiltonian now includes terms that 
reflect both the linear momentum of the particle in the potential field, and the spin angular 
momentum of the particle. The Schrodinger equation can be written as a two-component 
vector equation 
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where fL 



OJr, 



j m (J^ — s2 )- The subscripts + and — refer to the two levels of the atom 
(modelled in physics literature as the spin-up and spin-down states of a spin-1/2 particle). 
An applied field causes transitions between the eigenstates of the coupled spin-oscillator 
system. 

In order to understand the control of this system, we briefly present the a spin-0 particle 
in a quadratic potential with resonant frequency uj m being excited by a traveling wave of 
central frequency uj with controllable amplitude E(t). This system might be described by 
an equation of the form 

dip iu m ( d 2 x 



— x \ ip — iE(t)jji cos(kx — ujt)ip, 



(37) 



dt 2 \dx 2 

where \i is a dipole operator that describes the transition couplings between various eigen- 
states of the system, and k is the wavevector of the applied field. Let us introduce the 
differential operators 
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These are the familiar creation and annihilation operators of the quantized harmonic oscil- 
lator [35]. Introducing the differential operators corresponding the position of the particle 
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in the harmonic potential, and using the notation established in the previous subsection, 



Eq. 37 becomes 



dip 
~dt 



-iu m [ N + - ) ij) 



-i (^E(t)ficos(ut — kxo(a + at)j ip. 



(40) 



The product of k, the wavelength of the light, and xq, the amplitude of the zero-point motion 
of the particle in the harmonic potential (or the spatial extent of the ground state harmonic 
oscillator wave function) is the Lamb-Dicke parameter rj. Recall that the driving field is not 
resonant with the characteristic frequency of the harmonic potential. Therefore the above 
equation can be simplified only with the knowledge of the driving frequency, and the matrix 
elements of the operator /i that can be coupled by that driving frequency due to energy 
conservation. 



Returning to the system described by Eq. 36 , an applied field causes transitions between 
the eigenstates of the coupled spin-oscillator system. The amplitudes corresponding to 
the fields that cause the carrier and red transitions are dubbed E c and E r respectively. 
When both fields are applied simultaneously, ignoring overall additive constants, the control 
equations take on the form 
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H c E c (t) + fi r E r (t) 
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(41) 



where, fi Cjr is short for /xcos(cj C)r t — r}(a + al)). Here /i is the dipole operator that describe the 
strengths of the transitions between the two-level system eigenstates coupled by the fields. 
Examining the cosine terms in the off-diagonals, these can be expanded as 



cos 



(cut — r](a + a*)) = | (exp(zu;t) exp(— ir](a + at))) 
+ exp(— iut) exp(irj(a + at))) 



(42) 



Thinking in terms of the field adding energy to the ionic state (two level system), the term 
exp(— iut) adds energy fko to a state while the term exp(iut) removes energy hw from a 
state. For a true two-level system, one cannot add energy to the excited state, nor remove 
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energy from the ground state. So, keeping only the energy conserving terms (rotating wave 
approximation), the control equations become 
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(43) 



where /2 Cjr is short for /i exp(— io; c>r t + z^(a + a*)). 

In order to write the evolution equation in matrix form, the eigenstates are ordered as 



I T,o), | T,i>, .. 

block matrix form as: 



|,0), | .... The drift matrix of this system can be written in 



V 











=f- + Nu m 



(44) 



where, N is the previously defined number operator. The control matrices are not easily 
expressed in this basis. 

The evolution equation takes on a much simpler and elegant form in the interaction 
picture. Also, the transition matrices are described compactly by using the three Pauli 
matrices cr x , a y , and <j z . For example, the spin-flip transition is described by the operator 
°+ — °z + i&y and its hermitian conjugate cr_. In a coordinate system defined by the 
eigenbasis of the coupled system (\S,n), where S refers to the spin state and n to the 
harmonic oscillator number state), a general matrix element of the control Hamiltonian 
(S'n'\Hj\Sn) may be written as [TTf |4"3"]: 

'(S'\a+\S)(n'\exp(i(n(a + <J)))\n) 



-i2Re 



-i2Re 



(5V+|5)exp(-r ? /2),M ^t'^L^^) 



The symbol n> refers to the larger of n and n', and n < refers to the smaller of n and n'. 
Here L"(x) is the associated Laguerre polynomial. When the applied field of frequency uj c 
connects states | |,n) and | |, n) (carrier transitions), n' = n, and when the applied field of 
frequency uj r connects states | |,n) and | f,n — 1) (red sideband transitions), n' — n — 1. 
The matrix elements are zero for all other values of n! . 
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Derivation of the evolution equation for the trapped-electron qubit system 

We consider a trapped-electron qubit which can be well-modelled as a spin-half particle 
coupled to two harmonic oscillators. As detailed in Ref. |^6] . the relevant quantum states 
are the spin states of the electron (a spin-half particle), and the harmonic oscillator number 
states of the cyclotron and axial harmonic potentials. As before, the evolution is described 
by 



dtp+ 
dt 

dtp- 

at 



(45) 



u' c N c + u z N z + f 



(46) 



o vt tot V- 

where Vt tot = uj c (— + (3x^j +uj z (— -jj^ + xfj +^ s - Here, u s refers to the angular frequency 
associated with the spin-flip of the electron. The parameter (3 indicates that the cyclotron 
harmonic potential (indicated by subscript c) and the axial harmonic potential (indicated 
by subscript z) are of different strengths. 

In the basis of the harmonic oscillator eigenstates, the drift matrix can be written as 



oj' c N c + u z N z - f 

where the frequencies u c and u z refer to the frequencies of the cyclotron and axial harmonic 
oscillators, and uj s refers to the spin flip transition frequency. 

The energy eigenstates can be written as \nlj) where to refers to the number state of the 
cyclotron oscillator, I refers to the number state of the axial oscillator, and j refers to the 
spin state (up or down). The system is addressed by three control fields. The spin qubit is 
controlled using a small transverse magnetic field that connects states \ n I \) and \n I |). 
The spin-axial transition is controlled by a travelling magnetic field that connects states 
| n / |) and | to I + 1 j). A spin-cyclotron transition is controlled by an oscillating magnetic 
field with a spatial gradient that connects states \n I j) and \n — l I f ). The transfer graph 
in Figj5] clearly indicates that the three fields transitively connect all the eigenstates of the 
spin-axial-cyclotron system. 

As in the case of the trapped-ion, the evolution equation is clearer in the interaction 
picture. The evolution equation is given by 

z = (n s (t)B s + n sa (t)B sa + n sc (t)B sc )z. (47) 

Here Q s {t) is the control parameter for the spin flip transition, Q sz (t) is the control parameter 
for the spin axial transition, and fi sc (t) is the control parameter for the spin cyclotron 
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transition. The spin-flip transition operator B s = — i/2 (<7+ exp(— i<p) + er_ exp(i<f>)). The 
spin-axial transition operator 

B sa = -i(a+ exp(i(77(4e iw * t + a^e"^*) - A • t - 0)) 

(48) 

+ <r_ exp(-i(r7(4e^* + a z e" w ) - A • t - 0))) 

where the characteristic frequency of the transition is uo sa = uj s +uj z , the detuning A = uj—uj s , 
<p is an experimentally variable phase, and r\ is the Lamb-Dicke parameter. Here a z and 
a\ are annhilation and creation operators respectively associated with the axial harmonic 
oscillator. For uj = u sa (when the applied field is at resonance with the spin-axial transition) 
only the levels \il J.) and \n(l + 1) f ) are connected. The spin-cyclotron transition operator 
B sc = —i{a + a c e~ 1 ^ + o"_aJe*^) , where u sc = uo s — uu' c is the characteristic frequency associated 
with the spin cyclotron transition. Here a c and a\ are annhilation and creation operators 
respectively associated with the cyclotron harmonic oscillator. When the applied field is in 
resonance with the spin-cyclotron transition u = u sc , only the levels \(n — 1)1 |) are 

connected. 
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